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摘 要 : 激 波 计算 不 稳定 问题 是 激 波 


格式 ， 对 这 
的 多 维 性 ， 而 保留 了 
的 修正 方案 ， 改 进 了 


问题 进行 了 机 理 丰 


中 图 分 类 号 : 0354 


的 广泛 关注 与 看 


忆 在 了 


动量 插值 


机 秆 


网 表明 这 一 改进 上 


| 


上 修正 速度 中 的 压力 梯 
的 Roe 格式 可 以 有 效 的 消除 激 波 不 稳定 现象 。 
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F 激 波 计算 格式 未 能 真 ] 
度 项 。 通 过 这 一 机 型 
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有 获 格式 的 一 个 重要 缺陷 ， 多 年 来 受到 学 术 界 上 
究 ， 提 出 了 新 的 观点 ， 认 为 造成 这 一 问题 的 原 
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Abstract: The numerical shock instability is a famous and important problem for the shock-capturing 


scheme. For the classical Roe scheme, the mechanism of this problem is researched by the paper. A new 


viewpoint is proposed to explain the reason. In detail of the explanation, flow multi-dimensional property is 


not really reflected by the shock-capturing scheme, and thus a mechanism similar to the momentum 


interpolation method, the term of pressure gradient in the modified interface velocity, 1s held although it 


should disappear. Based on this viewpoint for the mechanism, the improved Roe scheme is proposed by a 


concise method. The numerical examples show that the improved Roe scheme can cure the shock instability. 
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0 引言 


近 几 十 年 来 ， 激 波 捕获 格式 
进展 ， 发 展 了 诸如 Roel、 
的 系列 格式 ， 基 本 解决 了 激 波 能 否 


难题 , 是 现代 计算 流 


本 力学 (CFD ) 辉煌 
然而 ， 对 于 计算 激 波 ， 以 上 激 波 


究 取 得 了 巨大 的 
SM 等 著名 
自动 捕获 计算 的 
的 成 就 之 一 。 
获 格式 仍然 


HLIL 户 、AU 


存在 一 些 缺 陷 未 能 解决 ， 如 红斑 、 马 赫 杆 、 奇 偶 失 
联 等 激 波 不 稳定 现象 外 。 对 于 这 些 激 波 不 稳定 现象 ， 


目 


(1) 格式 混合 方法 


前 的 解决 方式 主要 有 三 类 : 


研究 表明 高 精度 低 耗 散 的 激 波 捕获 格式 如 Roe 
格式 往往 受制 于 激 波 不 稳定 现象 ， 而 一 些 较 粗粮 自 
高 耗 散 格式 如 HLL 格式 则 没有 这 个 问题 。 因 此 一 种 
思路 是 将 这 两 类 格式 通过 一 个 压力 梯度 探测 函数 


WEAYy 


联 


域 使 
定 现象 的 机 型 


另 一 种 方法 将 激 波 计算 不 稳定 现象 归咎 于 
求解 器 (也 就 是 激 波 


日 


j 低 耗 散 格式 。 但 这 


四 


难 o 


黎 曼 求解 器 


方法 未 能 触及 激 波 不 稳 
EE， 探测 函数 的 构造 也 较为 
(2) 旋 车 


E 广 。 


是 
条 件 下 才 


一 维 的 简单 


i 获 格式 ) 在 多 维 情况 下 
也 就 是 说 ， 黎 曼 求解 器 是 


导出 来 的 ， 而 在 多 维 条 件 下 ， 通 常 的 


是 将 


性 。 因 


[5 


维 歼 曼 求解 器 简单 所 
此 ， 旋 车 
寻找 特征 方向 反映 流动 的 多 维 性 ， 如 旋转 Roe 格式 
， 从 而 可 以 消除 激 波 不 稳定 现象 。 


黎 显 
信安 


E 广 ， 未 真正 考虑 多 
求解 器 试图 改进 这 一 点 


党 要 


(3) 压力 梯度 修正 方法 


黎 曼 
仅仅 
一 维 
改 法 
维特 
通过 


第 三 种 方法 认为 激 波 计算 不 稳定 的 原因 在 于 格 
式 耗 散 中 的 压力 梯度 项 ， 从 而 通过 压力 梯度 探测 因 


子 进 行 修正 ， 在 压力 梯度 变化 大 的 地 方 取 消 压力 梯 


接 起 来 号 , 即 在 激 波 区 域 


使 用 高 耗 散 格式 而 在 


[x 


只 已 
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也 是 有 效 的 。 

本 文 提 出 了 一 种 对 激 波 计算 不 稳定 机 理 新 的 解 
释 ， 并 应 用 这 一 解释 针对 Roe 格式 给 出 了 简洁 有 效 
的 修正 方法 。 这 一 解释 与 方法 可 以 视 作 结 合 了 以 上 
所 述 第 二 类 与 第 三 类 方法 。 


1 主导 方程 及 格式 通用 形式 


Euler 方程 为 例 ， 流 体力 学 主导 方程 为 : 
60Q 9 oF ,ac =0 (1) 
Ot Ox Oy 
加 
2 
其 中 Q-|A | 为 守恒 变量 ，E -| 人 ?了 
pv puv 
pE u(pE+p) 
ey 
puv cr 
G=| ， 为 对 流通 量 。 
pvi+p 
[v(pE+p) 
对 于 大 多 数 数值 离散 格式 ， 对 流 数值 通 量 可 以 
用 如 下 通用 方式 表达 : 


™ 
妆 
mm 


(2) 


其 中 五 代表 中 心 差分 ，F, 代表 数值 粘性 ， 即 对 
通 量 由 中 心 差分 与 数值 耗 散 构 成 。 由 于 格式 
的 区 别 主要 在 数值 粘性 项 上 ， 本 文 将 只 讨论 EE 。 
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这 里 c 为 音速 ，U =nwu+nyv 为 界面 法 向 速度 ， 


n.、n, 为 界面 法 向 向 量 在 坐标 朋 
根据 文献 [7] 给 出 的 激 波 捕获 格式 统一 表达 方 


上 的 分 量 。 


式 ， 三 维 条 件 下 Roe 格式 的 耗 散 项 可 以 展开 为 如 下 
形式 : 
Ap [0 | [Ip 
oo | jp 
Fe = El A(pv) |+6pln, |+oU| pv (5) 
(pw) n, pw 
[A(pE) LV] Le® 
方程 (5) 右 端 项 具有 明确 的 物理 意义 ， 第 一 项 
为 基本 的 迎风 耗 散 项 ， 第 二 项 为 界面 修正 压力 ， 第 
三 项 为 界面 修正 速度 。 这 三 项 的 表达 式 如 下 : 
<=4= 人 O| (6) 
6p= 和 二 cp + [4 全 : [uap A(pU)| (7) 
中- p+ 入 名 -lvUAp-A( A(pU)] (8) 
p 2 
这 里 : 
1 vy 
A | 2 Ap wa(pu) walpr)* AlpE)| 
-人 
3 
(9) 


如 文献 [8] 所 述 , 第 三 项 5U 中 压力 梯度 项 的 机 制 


2 Roe 格式 及 其 展开 


类 似 动量 插值 ， 


在 低 马赫 数 条 件 下 起 到 抑制 压力 速 


度 失 侦 的 作用 。 


本 将 基于 


Roe 格式 进行 分 析 。 为 简便 起 见 


了 多 维 高 马赫 数 条 件 下 激 波 计算 不 稳定 现象 ， 将 在 


不 失 一 般 性 ， 本 文 的 讨论 只 针对 直角 正 交 网 格 i+ 


交接 面 。 
Roe 格式 的 数值 粘性 项 可 以 表达 为 如 下 形式 : 


FR --5R A 


di 


R (Qu —Q,) (3) 


| 
计 二 
2 


中 有 是 50 的 右 特征 短 阵 ，A 是 相应 的 特征 
值 对 角 和 矩阵 ， 其 中 特征 值 是 : 


= 上 和 = 区 -中 ， 


4 =|U +d| (4) 


下 节 讨 论 。 
3 激 波 计算 不 稳定 现象 机 理 分 析 与 改进 


如 前 所 述 ，6U 中 压力 梯度 项 的 机 制 类 似 动量 插 
直 ， 在 低 马 赫 数 计算 中 是 非常 必要 的 ， 否 则 将 导致 
压力 振荡 甚至 计算 发 散 。 然 而 ， 动 量 插值 方法 的 理 
论 核心 ， 在 于 界面 速度 中 的 压力 梯度 项 采用 计算 获 
得 ， 而 界面 速度 中 的 其 它 项 通过 插值 获得 。 也 就 是 
说 ， 本 该 由 插值 获得 的 压力 梯度 项 ， 在 动量 插值 法 
中 是 以 计算 获得 的 ， 因 此 与 其 它 项 并 不 匹配 。 

对 于 低 马赫 数 流动 ， 由 于 流动 的 非 线 性 不 强 ， 
因此 这 一 做 法 并 不 会 带 来 大 的 匹配 误差 ， 反 而 能 够 


~ 


抑制 棋盘 型 压力 的 发 生 。 然 而 ， 对 于 高 马赫 数 流动 ， 
一 方面 棋盘 型 压力 本 身 已 不 再 出 现 ， 另 一 方面 ， 由 
于 流动 的 强 非 线 性 ， 动 量 插值 法 造成 的 匹配 误差 已 
不 可 忽视 。 因 此 ， 在 高 马赫 数 流动 计算 时 ， 不 应 该 
再 有 动量 插值 法 或 其 类 似 的 机 制 。 也 就 是 说 ， 当 流 
动 马赫 数 趋 近 于 1 或 大 于 等 于 1 时 ，6U 不 应 该 再 包 
含 压力 梯度 项 ， 即 应 满足 以 下 条 件 : 


当 |M|>1 时 ，5U =0 or 5U = F(AU) (10) 
下 面 针对 Roe 格式 分 析 56U ， 即 公式 (8) 可 以 
进一步 展开 为 如 下 形式 : 
到 -| 全 4 名 -A (11) 
2 DC 2c 
可 以 看 到 ， 当 局 一 0 时 ， 公 式 〈11) 变 为 : 
i + 二 AU f (Ap, AU) (1 
DC 
se 公式 〈11) 变 为 : 
6U =AU= f(AU) (13) 
在 一 维 条 件 下 ， 可 以 看 到 ， 公 式 (12) 中 Roe 


格式 为 低 马赫 数 流动 提供 了 类 动量 插值 机 制 ， 而 公 
式 (13) 中 这 一 机 制 消失 ， 因 此 完全 符合 公式 (11)， 
不 存在 激 波 计算 不 稳定 的 问题 。 然 而 ， 在 多 维 条 件 
下 ， 存 在 这 样 一 种 可 能 性 ， 即 : 

U—0 但 M >1 (14) 

对 于 这 种 情况 ,从 马赫 数 的 角度 看 要 求 6U 满足 
公式 (10)， 但 此 时 6U 实际 上 是 由 公式 (12) 决定 
的 。 因 此 ， 本 文 认 为 造成 激 波 计算 不 稳定 现象 的 原 
因 ， 在 于 由 于 格式 未 考虑 流动 的 多 维 性 ， 而 保留 了 
本 应 消失 的 类 动量 插值 机 制 。 

从 这 一 观点 出 发 ， 修 正 Roe 格式 的 方法 也 将 非 
常 清楚 了 , 即将 具有 多 维 含义 的 马赫 数 引入 6U 中 的 
压力 梯度 项 ， 例 如 : 


Me 到 + 和 4 jp Dr] 


人 -和 [UAp-A(p0) 


(15) 


其 中 FM ) 为 马赫 数 的 函数 ， 构 造 方式 可 以 很 


多 ， 较 好 的 方式 可 以 参见 文献 [9]。 本 文 为 机 理 清晰 
起 见 ， 采 用 最 简单 的 定义 方式 如 下 : 


f(M)=min(M,!) (16) 


此 时 ,公式 (14) 所 引发 的 问题 即 可 获得 解决 。 


对 于 激 波 计算 不 稳定 问题 ， 本 文 所 提出 的 机 制 
解释 与 修正 方法 ， 与 引言 中 所 述 第 二 类 旋转 黎 曼 求 
解 器 方法 相 比 ， 同 样 认 为 是 由 于 流动 的 多 维 性 引起 
的 ， 但 具体 机 制 与 解决 方式 完全 不 同 ; 与 引言 中 所 
述 第 三 类 压力 梯度 修正 方法 相 比 ， 同 样 认为 原因 在 
于 压力 梯度 项 ， 但 指出 只 有 有 具有 类 动量 插值 性 质 的 
压力 梯度 项 才 是 问题 真正 的 根源 ， 并 且 修 正 函数 应 
基于 马赫 数 。 


4 数值 验证 


为 了 验证 修正 公式 (15) 的 有 效 性 ， 下 面 给 出 
经 典 的 三 个 算 例 。 
4. 1 运动 激 波 不 稳定 奇偶 失 联 〉 算 例 

该 算 例 中 ， 一 个 二 维 管道 中 ， 一 道 激 波 在 以 马 


赫 数 6 的 速度 向 前 运动 。 网 格 等 间距 为 1 并 且 平 行 
分 布 ， 但 管道 中 心 网 格 线 奇 数 点 与 偶数 点 坐标 分 别 
按 平行 线 坐标 上 移 与 下 移 0.0001. 

从 图 1 中 可 以 看 到 ， 使 用 经 典 Roe 格式 计算 该 
运动 激 波 ， 在 一 定 迭 代步 后 从 激 波 中 部 出 现 扰动 并 
逐步 放大 ， 最 终 无 法 保持 。 而 图 2 表明 ， 使 用 本 文 
改进 后 的 Roe 格式 , 即 采 用 公式 (15) 所 定义 的 6U ， 
可 以 很 好 的 保持 该 激 波 的 运动 。 


[ 


(a) 6000 和 迭代 步 


人 


(b) 10000 和 迭代 步 


Cc) 14000 迭代 步 
图 1 经 典 Roe 格式 计算 结果 
Fig.1 Results by the classical Roe Scheme 
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(a) 6000 和 迭代 步 


| 


(b) 10000 和 迭代 步 


上 


(c) 14000 帮 代 步 

图 2 本文 改进 Roe 格式 计算 结果 

Fig.2 Results by the improved Roe Scheme 
4. 2 双 马 赫 数 反射 问题 

一 道 运 动 马赫 数 为 10 的 激 波 , 以 60” 和 角 扫 过 平 
板 ， 形成 复杂 的 激 波 系 。 图 3 与 图 4 给 出 了 0.2 秒 时 
计算 所 得 的 密度 云图 。 可 以 看 到 ， 采 用 经 典 Roe 格 
式 时 ， 激 波 系 中 本 应 与 平板 近乎 垂直 的 激 波 变 得 倾 
斜 并 且 出 现 了 非 物理 的 三 贫 点 现象 ， 即 马赫 杆 问题 。 
这 也 是 激 波 计算 不 稳定 的 经 典 现 象 之 一 。 而 本 文 改 
进 的 Roe 格式 则 较 好 的 避免 了 这 一 问题 。 


图 3 经 典 Roe 格式 
Fig.3 The classical Roe scheme 


CC 
图 4 本文 改进 Roe 格式 


Fig.4 The improved Roe Scheme 
4.3 高 超 音 钝 体 绕 流 问题 


图 5 经 典 Roe 格 式 ”图 6 本 文 改进 Roe 格 式 
Fig.3 The classical Roe Fig.6 The improved Roe 
马赫 数 为 20 的 高 超 音 来 流 绕 流 通过 如 图 5 与 图 

6 所 示 的 圆 球 钝 体 ， 并 在 球体 前 缘 形成 了 强 脱 体 激 

波 。 图 5 中 采用 经 典 Roe 格式 所 获得 的 结果 ， 在 正 

前 方位 置 的 激 波 出 现 扭曲 现象 ， 即 激 波 计算 不 稳定 


的 经 典 现 象 之 一 红斑 问题 。 而 图 6 中 采用 本 文 改 进 


Roe 格式 的 结果 则 校正 了 这 一 问题 。 


5 结论 


提 昌 


本 文 对 激 波 计算 不 稳定 现象 进行 了 机 理 分 析 ， 


4 了 新 的 看 法 ， 即 认为 原因 在 于 由 于 激 波 计算 格 


式 未 真正 考虑 流动 的 多 维 性 ， 而 保留 了 本 应 消失 的 
类 动量 插值 机 制 。 通 过 这 一 认识 ， 提 出 了 简明 的 修 


正方 案 。 数 值 算 例 表 明了 这 一 修正 方法 的 有 效 性 。 
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